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We systematically study the question of identification and consistent inclusion of the radion, 
within the Lagrangian approach, in a two brane Randall-Sundrum model. Exploiting the symmetry 
properties of the theory, we show how the radion can be identified unambiguously and give the action 
to all orders in the radion field and the metric. Using the background field method, we expand the 
theory to quadratic orders in the fields. We show that the most general classical solutions, for 
the induced metric on the branes in the case of a constant radion and a factorizable 4D metric, 
correspond to Einstein spaces. We discuss extensively the diagonalization of the quadratic action. 
Furthermore, we obtain the 4-dimensional effective theory from this and study the question of the 
spectrum as well as the couplings in these theories. 



I. INTRODUCTION 

The idea that we may live in a world of more than four space-time dimensions is certainly not new , but in recent 
years it has received a lot of attention because of the observation that the extra compact dimensions need not all be 
of the order of the Planck length scale. A large extra dimension of the order of a (TeV)^ 1 || provides an alternate 
solution to the hierarchy problem (the desert scenario) , namely, it explains in a natural manner the large ratio for the 
Planck scale to the electro-weak scale. There are basically two interesting approaches to solving this problem. The 
first, due to Arkani-Hamed, Dimopoulos and Dvali is based on the standard Kaluza-Klein approach, where the 
higher dimensional space is a direct product of the 4 dimensional space-time and a compact manifold. The second, 
due to Randall and Sundrum (RS) P,H, uses an alternate and interesting scenario, where the higher dimensional 
space-time is not factorizable. 

The simplest RS model can be envisaged as follows. Let us consider a five dimensional space-time manifold, with 
the extra coordinate taking values on an orbifold S 1 /Z2, with two fixed points, which are chosen to be at, say, and 
ir. At each of the two fixed points of the orbifold, there is a singular source - a 3-brane - which carries a constant, 
nonzero tension. In addition, a negative five dimensional cosmological constant is also assumed to be present in the 
bulk, which leads to a five-dimensional metric that is a slice of the AdS space || . Since the warp factor turns out to 
be an exponential function depending on the extra coordinate, the ratio of the values of the warp factor at the two 
branes, for such a metric, naturally leads to the ratio of the Planck to TeV scales without the need for an extremely 
large extra dimension. We can think of our physical world as described by one of the two branes, conventionally 
chosen to be the one at ir. In the original RS model, the values of the five dimensional cosmological term and the 
tensions of the two branes are chosen in such a way, that the effective four dimensional cosmological term is zero and 
the induced metric is just the Minkowski metric. In this model, the physical distance between the two branes is an 
arbitrary constant r c , whose value cannot be determined from the classical analysis. The stabilization of the size of 
the extra dimension requires additional mechanisms j5|. Even when the five dimensional cosmological constant and 
the brane tensions are not fine tuned to cancel (as is the case in the original RS model), it was shown in Q that the 
model admits solutions with induced metrics of dS or AdS type. In such a case, the distance between the two branes 
can be determined in terms of the parameters of the model. 

The classical solutions of the RS model, as discussed above, correspond only to the vacuum sector of the theory 
when there is no matter present. One can, of course, add matter to the theory, which is assumed to be located on the 
physical brane. In principle, one can also add bulk matter. In addition, one can also analyze the fluctuations of the 
five dimensional gravity around the vacuum solution, as well as the fluctuation of the distance between the branes. 
In this paper, we will not consider the inclusion of matter, and will investigate only the contributions arising from 
the fluctuations of the metric and the geometry. To fix the terminology, let us note that the four dimensional scalar 
field, which describes the distance between the two branes, is called the radion. The vacuum expectation value of this 
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field is expected to determine the constant r c , although we will not worry about the mechanism that leads to this 
vacuum expectation value. The radion was first introduced in ||,^| , where fluctuations of gravity and the radion were 
simply added to the vacuum values of the metric and r c respectively. Although this is quite natural and is standard 
in quantum field theory, it was later pointed out in that the ansatz of |3|,|| does not solve the linearized Einstein 
equations. Instead, the authors of Q used an alternate form of the metric, that is correct to linear order in the radion 
fluctuations. 

The inclusion of the radion field was of great importance in the construction of the supersymmetric version of the 
RS model (|J^1, where the bosonic sector of the theory, in addition to other fields, contains the radion and the graviton 
(but not its KK tower) to all orders in these fields. The ansatz used in || coincides to linear order with that of 0, 
while that in corresponds to pM]. The two papers contain similar results, which are valid up to second order in 
the space-time derivatives. The inclusion of the radion in the case of dS and AdS branes has also been discussed in 
p0| , using an ansatz of the form of [pLEJ . 

We believe that, at present, we do not have a complete understanding of the radion in the RS model and it is the 
proper identification and the consistent inclusion of the radion field in an arbitrary two brane model which is the 
purpose of this paper. We derive the complete five dimensional Lagrangian for the RS model, including the radion and 
the graviton (with its KK tower), and study systematically the question of the spectrum of the theory by restricting 
the Lagrangian up to second order in these fields. We follow the geometrical approach of Jjj, where the authors 
identified the radion field from an analysis of the symmetries of the Einstein equations in the bulk, and derived the 
form of the five dimensional metric up to linear order in the radion. The work of JtJ emphasizes the importance of 
the junction conditions which, as we will see, are of utmost importance in our investigation as well and are quite 
relevant for understanding the mixing between the radion and the graviton and its KK tower. However, in our opinion, 
there are some features in the analysis of that remain obscure and we study this problem systematically from the 
conventional Lagrangian approach. We formulate the RS model as the model of five dimensional space-time with 
two 3-branes at the fixed points of an orbifold, whose action is manifestly invariant under arbitrary five dimensional 
transformations, and we carefully elucidate the symmetry properties of the model. Such an approach makes the 
identification of the radion field unambiguous and free of assumptions. Subsequently, employing the background field 
method, we are able to make a very general and systematic investigation of the possible vacuum solutions in the RS 
model as well as give a satisfactory 51? Lagrangian description (that holds for any vacuum solution) of the radion and 
the graviton and its KK tower, up to second order in the fields, without mixing between them. 

The paper is organized as follows. In section 2, we describe the model, both in an interval as well as on the 
orbifold. Taking advantage of the symmetries of the theory, we parameterize the metric in a manner that makes 
the identification of the radion natural and is convenient for our subsequent discussions. The Lagrangian density 
of the theory is then expressed completely in terms of this parameterization to all orders in the field variables. In 
section 3, we use the background field method to expand the action as well as the junction conditions to quadratic 
order in the field variables. The classical equations are solved in a unified manner and we show that a factorizable 
background metric, in general, defines an Einstein space multiplied by a warp factor. We then discuss the question of 
diagonalization of the quadratic action as well as the boundary conditions in a systematic manner. In section 4, we 
rewrite this action as an effective action in 4-dimensions and discuss further the properties of this diagonalized theory. 
In section 5, we present a brief summary of our results. In appendix A, we describe some relations that are useful in 
the background field expansion of the theory while in appendix B, we discuss the Kaluza- Klein decomposition of the 
metric as well as some of the properties of the basis functions. 

II. THE MODEL 

Let us consider a 5-dimensional space-time manifold with signature (—,+,+,+,+), which is parameterized by 
the usual 4-dimensional space-time coordinates as well as a fifth coordinate that is bounded by the locations of two 
4-dimensional hypersurfaces. The hypersurfaces can be specified by equations of the forms 

f L (x,z) = 0, f R (x,z)=0 (1) 

The functions Jl,r are, in general, arbitrary except for the condition that everywhere on the boundary hypersurfaces, 
the normal fields, given by n = dfL,n(x, z)|/=o, ar e space-like. The locations of the two hypersurfaces can be 
determined by inverting ([!]), namely, 

zl = 4>l(x) , z R = 4>r{x) (2) 
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where 4>l,r(x) are not necessarily small. Therefore, z £ I z := [<Pl(x), 4>r(x)]. From (||) it is clear that under arbitrary 
AD transformations, z' — z , x' = x'(x), the functions 4>l,r transform as AD scalars. 

On this 5-D space-time manifold , one can define a theory which is invariant under arbitrary 5D transformations. 
The action has the form 

S = Stmik + Sl + Sr (3) 

where 

S bulk oc / d 5 V^G (-R^ — A + . . (4) 

./bulk V ' 

Si oc f d A Vyf^{-Vi + ...), i = L,R (5) 

Jfi=a 

Here, Gmn and G are the 5D metric and its determinant respectively, while g l mn is the AD induced metric on the 
respective boundary hypersurface. In this paper, we use the convention that capital roman letters denote 5D indices, 
while lower case roman letters represent AD indices. We also use the notations and definitions of throughout. A is 
a 5D cosmological constant, Vi are the tensions on each of the branes and the dots denote possible terms representing 
5-D as well as AD matter, which we ignore in the present study. 

To obtain the equations of motion, we need to vary the action (|^). As long as the positions of the boundaries are 
not specified (i.e. Jl,r are arbitrary), it is not enough to consider the variation of the metric alone. We must also 
extremize the action with respect to the "volume" where the theory is defined. To do so, it will be more convenient 
to make use of the 5D coordinate invariance of the model and choose a special coordinate system, where the volume 
does not depend on dynamical quantities such as fz,R and, consequently, is not subject to variations. 

Such a transformation can be done in two steps. First, we adopt Gaussian Normal (GN) coordinates with respect 
to the brane </>£. In this case, we have 4>l — constant, which we choose to be zero, namely, <pi, = 0. In this coordinate 
system, the second boundary is located at 4>r = 4>{x) and the 5D metric takes the form, G r „5 = , G55 = 1. Clearly, 
the function (f>(x) represents the physical distance between the two boundaries, and can be related to the radion field. 
Next, we rescale the fifth coordinate z = 4>(x)t/ir , t £ I t := [0, it] (t does not correspond to time which is denoted by 
x ). Introducing a new function 4>(x) := 4>(x)/ir, we have the two boundaries at the fixed end points t = and t = ir 
respectively and the 51? volume takes the form: M 4 x I t . (A more general transformation of this type has already 
been discussed in Jl2|] .) Note that neither the above transformations nor the field <fi(x), which we call radion, have to 
be small. It is worth noting here that <f>(x) is a four dimensional scalar field as the radion should be and it has no 
dependence on the extra coordinate, a reflection of the fact that there is no Kaluza-Klein tower for the radion. 

In these coordinates, the 5D metric takes the form 

(Qran 
(6) 
N m 4> 2 (x) ) 

and the inverse metric is determined to be 

G MN = (7) 
V w iv / 

where the raising and the lowering of the AD indices, to, n, is done with the metric g mn (x,t) satisfying 
g m n(x,t)g nk (x,t) = 6^ n . Furthermore, N m (x,t) = |(0 2 (x)) !m and N(x,t) = <f> 2 (x) - N m (x,t)N m (x,t) where a 
comma denotes a derivative. We see that the metric in these coordinates has a form similar to that of Arnowitt- 
Deser-Misner |l3"| . Let us note that, in these coordinates, the induced metric on each of the boundaries is given 
by 

SmnfX) = 9mn(x,t)\ t = ti - (8) 

Furthermore, it is easy to show that 

det(G MN {x,t)) = N det (g mn (x, t)) . (9) 
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This defines the theory on an interval t G If = [0, w]. Let us next replace the interval, It, by the orbifold S 1 /Z^- In 
order to do that, we extend the variable < t < tt to the interval [— 7T, 7t], imposing the additional symmetry t — > — i. 
As a result, we need to replace t by r = |f| in the metric @ so that g mn (x,T) is a symmetric and non-degenerate, 
but otherwise arbitrary tensor and, 



N m {x,T) = -(4> 2 (z)), m 
A(x,r) = 4> 2 (x)-N m (x,T)N m (x,T) 



(10) 



Throughout the paper, we will denote by a dot (e.g. a) differentiation with respect to r, while a prime (e.g. a') 
denotes differentiation with respect to t. We denote the AD partial derivatives by a comma (e.g. a, m ), while a 
covariant derivative is represented by a semicolon. 
The action (0), on the orbifold, can be written as 



S — Sbuik + ^ Si 



i—0,7T 

d 4 x / dt V^G f -2M 3 i? (5) - A 



Sbulk ~ 

S t = I d*x I dt y/-g* (-V<) *(t - fc), 



= 0, 7T 



(11) 



where i?^ 5 - 1 is the Ricci scalar constructed from the 5D metric Gmn given in (^3j) and generalized to the orbifold, 
as discussed above. In studying the solutions of this theory, the standard approach would be to solve the Einstein 
equations following from the action ([n]). However, we note that, since the 5D metric (JgJ) is already decomposed into 
a 4 + 1 form because of the special coordinate choice we have made, it is natural to recast the action first in terms 
of the AD metric, g m n(x,T), and the radion field, 4>(x), using this parameterization. After some lengthy algebra, we 
obtain 



5* 



-2M' 



d A x 



i— 0.7T 



(12) 



where 



A 



N m N r, 
N 



^ [g' mn (g mn y + [g mn g' mn ) 2 ]) + c Mgh 



c 



N 



high 



v2 {g mn N°N s , n {4> 2 ), m - A%A m (^ 2 ),, 
+ 2N m N n [N k . k N m;n - N k . m N k . n ]) 



Vi 
2AP 



S(t-U) 



(13) 

(14) 
(15) 



In the above expression, Rmli is the Ricci tensor constructed from the metric g mn , and a semicolon denotes the 
covariant derivative with respect to this metric. Because of the dependence on the fifth coordinate, all these operations 
are to be carried out for fixed values of t. The tensor nature of each of the terms, with respect to AD coordinate 
transformations, can be read off easily. 

Here, we would like to emphasize the fact that, although our starting theory was manifestly invariant under 5D 
coordinate transformations, only a residual AD symmetry survives in ( |l2|) because of the special choice of coordinates 
made. Thus, action ( |l2| ) is invariant under arbitrary 4£)-coordinate reparameterizations, x — > x'(x). Let us also note 
that, since all components of the metric, Gmn, are even functions of the extra coordinate, it follows from the chain 
rule, F'(t) = F(t)t' , that expressions linear in ^ derivative are odd and, therefore, do not contribute to the action. 
We have omitted such expressions in (13). Furthermore, we have also dropped surface terms arising from integration 
by parts. This, therefore, gives the complete action involving the graviton (and its KK tower) and the radion. 
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III. BACKGROUND METHOD 



To study the solutions of this theory systematically, we use the background field method [|L4 15 (see in particular 
section 4 in [|l5|). We write each of the field variables as a sum of a classical (possibly large) background and a small 
fluctuation, and expand the action up to quadratic terms in the fluctuations. For the metric we write 

g mn ( X ,T) = g mn (x,T) - h" m (x,T) + h m k (x,T)h kn (x,T) (16) 

while, for the radion, we use the decomposition 

4>{x) = r(x) + r{x)d{x) (17) 

where g m n{x,T) and r{x) are the background fields which are assumed to satisfy the classical field equations, while 
h mn and d are the corresponding fluctuations (note that we have introduced d as a dimensionless field). The new 
(background) metric satisfies g m kg kn = ^mt an d from now on all indices are lowered and raised with the background 
metric g mn and its inverse. For simplicity, we will refer to this as the metric. 

Our next task is to decompose the action up to second order in h and d. The zeroth order action leads to the 
classical equations for the background fields. The linear terms vanish, since the backgrounds satisfy the classical 
equations and, therefore, the first nontrivial term in the expansion of the action corresponds to the quadratic part of 
the action in the h and the d fields. We note here that our interest lies in the case where the classical solutions lead 
to r(x) = r c = constant primarily for two reasons. Such a case would correspond to have the highest symmetry and 
the calculations will be much simpler. When r(x) is a constant, it is easy to see that the terms in (|lj) would contain 
terms that are at least third order in d(x) and, therefore, will not be relevant for our analysis. It is for this reason 
that we have separated out these terms in (|l^) and we will neglect this term in the rest of our discussions. On the 
other hand, if non-constant solutions for r(x) are of interest, then, the term in (|lj) will contribute to the expansion 
up to quadratic order. 

Before carrying out the expansion, however, let us discuss the Israel junction conditions fl6[| that are crucial in 
analyzing solutions in this theory. Let us note that we are considering a theory with singular sources and, in such a 
case, boundary conditions on the branes are extremely important in determining solutions. The junction conditions 
for the metric g mn , can be derived following jlTj , and in the present case of an orbifold read as 

V, 



g mn (x,T)\ t= t t = -£(t) ^ 3 \J N(x, rjg mn (x, T)\ t =u (18) 



where 



(19) 

Decomposing the metric and the radion as in (|l6|), ( fi"7j ) and using ( |l0| ) as well as the results from appendix A, we get 
order by order in the fluctuations 

9mn( x > T )\t=ti = -SWj^j^SmnfoT)!*^ (20) 
h mn {z,T)\t=ti = -t(t)j^j^(h m n(x,T)+d(x)g mn {x,T))\ t=ti (21) 
and so on. The equations following from the action ( |l2| ) need to be solved subject to these boundary conditions. 



A. Classical Action 

The zeroth order action (which does not contain d or h) is obtained from ( |l2| ) with g and cf> replaced by g and r 
respectively. The equations of motion for the two fields can be obtained by varying with respect to g and r which 
leads to two equations - one tensor and one scalar. Let us note parenthetically that <5G m 5 (see (^)) is not independent 
and, consequently, the m5 Einstein equation would not be independent. 

The zeroth order action or the classical action is easily obtained from (fl^) to have the form 
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5 (o) 



-2M 3 / d 4 x 



dty~9 r R + r 



A 



2Af 3 



- 5: [<4 n (<T")' + GT-SL.) 3 ] + E " **) ) 

4=0, 7T / 



The Euler-Lagrange equations, following from this action, have the form (upon setting r(x) = r c ) 

= R 



X +^[9' mn (9 mn Y + (9 mn 9 , mn ) 2 ] 







R 



2M 3 4r? 
1 



A 



ink 



2 R ~ 9mk 4AP 



1 



9mk 



4g'^ k - Ag ab 9ab9mk 



-^ ab 9' am 9bk ~ H~9 ab y~g'a b ~g m k + 2~ g al % b g' mk - (g ab g' ab ) g mk 



E 



Ar r M 3 



S(t — ti)g m k 



(22) 



(23) 



(24) 



As expected, equation ( p4| ) contains singular terms proportional to (^-functions. However, let us note that these 
equations, which are written in terms of derivatives with respect to t € [— tt, tt], can be simplified if we express them 
in terms of the variable r £ [0, 7r] (Recall that all the field variables depend on r.). We note that 



3mn( r ) = 9n 



gjnn (y) 9mn 9r> 



Since r' 2 = 1 , in the classical equations, we can simply replace terms quadratic in single t derivatives (primes) by those 
with r derivatives (dots). However, in terms with double t derivatives, the change to r derivatives would introduce 
new delta function singularities because of the r" term. On the other hand, we note from ( pfj| ) that g mn , on both 
the branes, is proportional to g mn so that these new singular terms cancel precisely the singular delta function terms 
already present in the equation (^). In fact, this is a general result, namely, singular terms arising from changing 
"primes" to "dots" exactly cancel the singular boundary terms that appear explicitly in the equations of motion, 
because of the junction conditions. We would, however, like to stress that, although the dot derivatives appear to 
be the natural ones in the equations of motion, it is more useful to have the prime derivatives in the action to avoid 
subtleties in integration by parts. 

Rewritten in terms of the r (dot) derivatives, the equations of motion have the form 



= R 



R 



A 1 

2M 3 + irf 
1 



g mn g + (g 



ran) 



rak 



-R g-mk 



A 



1 



rffmfe 



[ 4 5mfe - ^gabgrnk 



b ta m gbk - 3 <T SabSmfe + 2g ab g ab g mk - (g ab g ab ) g mk 



(25) 



(26) 



without any singular terms involving (^-functions. These are highly nonlinear equations which are clearly nontrivial 
to solve in general. Therefore, we look for a solution of the metric, <? mn , in the factorizable form, 



g m k(x,r) = a(r) g mk (x) 



Requiring g mk {x) g kn (x) = <5™ , we have 

1 



g mk 



- g 

a — 



Rmk — R m k. 

Using ( |27| ) and (|28|), equations ([25j) and ( |26| ) respectively take the forms 

2 



R=- R 

a 



= R + a 

= Rmk ~ 



A 

2M 3 

1 



4 ° 

A 



— R a , H — a , 

2 — —mk zj. —mk 



(27) 
(28) 

(29) 
(30) 



G 



where we have defined 



A 6 a 



M 3 



r? a 



(31) 



Contracting ( pOj ) with g mk , we obtain 

R = X (32) 



The left-hand-side of (|32|) is a function of x only, while the right hand side depends only on r. It follows, therefore, 
that A = constant = R, which, in turn, implies that the metric g mk defines an Einstein space 

&nk = \l mk (33) 

with constant scalar curvature A. It follows now from (|J) and (Pq) that the induced metric on each of the branes will 
define an Einstein space as well, with curvature R(t = ti) = a ^ t .y We also note that every solution, g mk , of ( |3"3| ) 
corresponding to a given A, will have the same warp factor a determined from (^9|) and (|3l|). 

The solutions for the warp factor a are already known. The flat case, g — r/ m k, was analyzed in ^j, while the 
dS and the AdS solutions were derived in fijl. It is clear from the preceding discussion that these solutions exhaust 
all possible cases for the warp factor and here we will briefly describe an alternate and unified derivation of these 
solutions for completeness. 

The equations (E9L) and (|3l|) determining the warp factor can be rewritten in a more convenient form 



£V = r " ( — i _ 

a) 3 \2AP a 



(34) 



Ar 2 

da — (a)" = a (35) 

We note from the above equations that a constant solution for a(r) exists only for the physically uninteresting case 
when A = = Vi, which we will not consider. When a ^ constant, equation (|35|) follows from ([m|), and, therefore, 
this is the only equation that we need to analyze. Furthermore, this equation must be supplemented by the boundary 
conditions (see (^0|)) 

a\ t=u = (36) 

The general solution of (|34|) is of trigonometric type for A > 0, of exponential type for A < and of polynomial 
type for A = 0. Keeping in mind the spirit of the RS model, where the aim is to solve the hierarchy problem without 
introducing large numbers, we will consider only the case A < here. Following H, we define 



A 



24M 3 

and parameterize the tensions on the branes, for convenience, as 

V , V. 



(37) 



^, -riH = -P* (38) 



24M 3 ' 24M 3 

where a, (3 are arbitrary parameters. In terms of these parameters, the equation for a reads 



Ar? 



= 4fc 2 r c 2 - ^ (39) 
a J 6a 



which needs to be solved subject to the boundary conditions 



a (r = 0) = -2akr Cl -(t = tt) = -20kr c (40) 



a a 
The solutions of (|39|) subject to (EJ) are easily determined to be 
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y^p- cosh 2 (fcr c (p — r)) when A > 0, 
a(r) = ^ exp(— 2fcr c r)) when A = 0, 
Y^ja sinh 2 (fcr c (p — r)) when A < 0, 

When A 7^ 0, r c is determined completely in terms of a and /3 as 

a-/3 



kr r 



iarcth 1 

7T VI 



a/3 



and the parameter p is defined through 



2kr 



2fcr, 



1 + a 
1 - a 
a + 1 



a 



1 



when A > 0, 
when A < 0, 



(41) 

(42) 

(43) 
(44) 



The range of the parameters, a and /3, can be obtained from three consistency relations. We note that r c > 0, which 
through ( [i"2"| ) implies that < ^E§» < 1- Furthermore, the right-hand-side of ( (43| ) must be positive, so that \a\ < 1 
for A > and \a\ > 1 for A < 0. Finally, we have a(r) > for any r and it follows from (|4l]) that p ^ [0, n] when 
A < 0. It is easy to show that, for the solutions in ([4l|), the right-hand-side of (39) is non-negative for any r £ [0, n] 
and, consequently, this does not introduce any further restriction on the parameters a and (3. The limiting case A = 
cannot be obtained from these and a solution exists only if the two tensions satisfy a = j3 = 1, as can be seen from 
(O) and (]4(]). In this case, the value of r c remains undetermined. 



B. Quadratic Action 

If the background fields satisfy the classical equations, then, the linear order terms in the expansion of the action 
vanish, which can be explicitly checked. The leading nontrivial correction to the action comes at the quadratic order. 
After some algebra, the part of the action (12), quadratic in h mn and d, can be determined to be (we give some useful 
relations for this expansion in appendix A) 

S^l = -2M 3 r c J d^x^lj dt (L(d) + L{dh) + L(h)) (45) 

where 



A 9 9 . ™ t 3 

(a 

aX 3 aa' 



L{d) = '-7r 2 c T 2 d> m d, m -^(a') 2 d 2 (46) 
4 ri 



L(dh) = a(h mk k -h' m )d m -^dh+^dti (47) 
V / ' 4 2 r~ 

L(k) = \ (±h ak . m h ak -' m - h a ^X m * + h m %h m - \k m h n )j (48) 

+£i (h! mk (h mk y - 0O 2 ) + ~ (2hf mk h mk - h'h) (49) 



° A ' 3(a ' n ^--^ (50) 



T 

4 rl 

-T, a ^(^ k hr k -h 2 )s(t- ti) (5i) 

i— 0,7T 

Keeping the factorization of the background metric in mind, we have introduced a new field h mk (x, r) through 

h mk = ah mk ; h mk =ah mk ; h = h mk g mk (52) 
All indices are raised and lowered with the background metric q , , which also defines the covariant derivatives. As 

. —mk' 

a result, let us note that 4D covariant derivatives simply commute with derivatives with respect to t (or r). 
The boundary condition satisfied by h mk is determined from (Ell) to be 
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h mk {x,T)\ t=u = -£{t)j^^ d g mk \ t=u (53) 

The most important feature to note from ( [45] ) is the mixing between the fields h mk and d. In order to understand the 
physical spectrum of the theory, we must, of course, diagonalize the second order Lagrangian density. Unfortunately, 
this is not the complete story since there is mixing between the two fields arising from the boundary conditions ( j53|) as 
well. Therefore, we need to diagonalize simultaneously the quadratic Lagrangian density, (|45|), as well as the junction 



conditions to linear order, (53). We note that higher order mixing, both in the action as well as in the junction 
conditions, would simply correspond to higher order interactions and is not of interest for our analysis. 

The presence of mixing, in the quadratic action, simply corresponds to the fact that h mk and d are not the 
appropriate field variables and the physical fields will be, in general, a linear combination of the two. However, we 
note that, unlike h mkl the field d, in some sense, carries a direct physical meaning - it labels the deviation of the 
distance between the two branes from the classical value r c . Consequently, we prefer to treat the field d as physical 
and look for a general redefinition of the metric fluctuation. The most general linear redefinition of h mn can be seen 
to be 

h mk (x,r) = X mk (x.T) + {\na + f{r))d{x) g_ mk + S{r)d, mk + Z( T )d%{x) g_ mk (54) 

where compatibility with (|5^) requires (see also (j36|)) 

X mk (x,r)\ t=u = (55) 
f(r)\ t=u = S(r)\ t=u = Z(r)\ t=u (56) 

Note that the term L(h) is form-invariant under the transformation (jHJ). However, L(d) and L(dh) do transform 
and it is easy to see that if Z ^ 0, such a transformation will generate mixing between X mk and d with fourth order 
derivatives. Since this does not help in diagonalization, we set Z = 0. Furthermore, we find that there is no such / 
and S obeying (|5^) which set the mixing terms, in the Lagrangian density, to zero. 

There is another possible solution to the mixing problem. We note here that it is the action that we need to 
diagonalize (subject to the boundary conditions). Therefore, if we know the KK decomposition for the field X mk , 
then the action ( |45| ) would be a 5D action for AD fields that depend only on the coordinate x, with t-dependent 
coefficients. It would, therefore, be sufficient to arrange the coefficient functions in front of the mixing terms to vanish 
when integrated over the extra dimension so that such terms can be removed from the action. 

To find a set of basis functions, in terms of which we can decompose the field X mk , we proceed as follows. Let 
us assume that the change ([34]) leads to a diagonalization of the action (it also diagonalizes the boundary conditions 
because of (pq)). Since the transformation (|I|) is linear and leaves the part of the action depending only on h mk 
form- invariant, the other two terms in the action can only lead to terms that are quadratic in the field d. As a result, 
the KK decomposition of X mk can be easily studied by analyzing the h mk part of the action and we can simply set 
d = in (|45|) and (^||) for this purpose. Having the basis functions determined in this way, X mk can be expanded in 
this basis (d(x) does not have a KK tower) in the action fl45|), and the resulting action diagonalized at each level of 
the KK tower. As we will see, this procedure works quite nicely. We discuss the actual solutions of ( (45| ) and (|53|) for 
d = in appendix B. 

The most convenient way to implement the program discussed above, is to first make the change of variables 

kmk{ x > T ) = x rnk(x,r) + (\na + f{r))d{x) g mk (57) 

with 

X mk (x,r)\ t=u = f( T )\ t=u = (58) 

The function /(r) is otherwise arbitrary. Since the new field X mk satisfies (p^), we can decompose it in the system 
of basis functions <I> Q (whose derivative vanishes at the boundaries, see appendix B for details) as 

X mk (x lT ) = (Q) A mfe (x)$«(r) (59) 

OL 

The action ( [45] ) can now be written in the form 

S (2) = -2M 3 r c f d 4 Xyf=g~ f dt (L(d) + L(dX) + L(X)) (60) 
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where 



L(d) = K(t) d< m d rn + 



M(t) 



(X m % - X' m ) d, m - -dX 



L(dX) = o(l + Ino + /) 
L(X) = fl (*«) 2 L PF {^X mk -ml) 



3a 2 
2^ 



f + 2-f dX 
a 



Here we have defined K and M as 



A-(r) = ^r c 2 r 2 -^a(ln a + /)(ln a + / + 2) 
M(r) = -6aa(lna + /)(lna + / + l) - 3o/(o/ + 2d(lno + /)) 
+ E a 2 ^(lna + /) 2 ^-t,) 



i=0,7T 



and represents the Pauli-Fierz term that we will discuss later (it is also discussed in appendix B). 
Because of (pq), we can also decompose / in terms of <I> Q as 



fir) = E ^» 



(61) 
(62) 
(63) 

(64) 
(65) 



(66) 



Using (BIO) and (B12), we can rewrite the mixing terms, (p2|), in the action in the form 
L(dX) = a(l + lna + p)$" ' " "" ' ' 



+ ^]a($ Q ) 2 (f +/ a )[( (Q) I m f fc - Wl' m )d, ro -^Wl 



^m 2 d (°>I 



where we have introduced the notations 



P = fo®° = constant, I a = / dt a In a $ Q , a > 1 

J — 7T 



(67) 



(68) 



and have used the equality J* dt alna$ Q = dt a($ Q ) 2 J a . 

It is impossible to find a set of constants, f a , which will make the mixing terms in the action vanish. As a result, 
we are forced to make a second transformation involving the second derivative of d, which is done at each level of the 
KK tower 



{a) X mk {x) = {a) H mk (x) + a a d ;mk , a>0 



(69) 



The numbers a a can be thought of as the coefficients in the decomposition of the function S(t), in (gl), in the basis 
functions namely, 



(70) 



Under the transformation (J69J), the Pauli-Fierz Lagrangian (see ( B18 )) becomes 

L PF {^X,m 2 a ) = L PF (^H,ml) - m 2 a a a ( ^H m ^ k - < a >ff' m ) d, 

I v.m \ 2 , ^ a a m a j,m 1 



(71) 



and the mixing term (|67|) takes the form 
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( 



(0) H mk _ (0) H ,m) d 



L(dX) = a(l + lna + p)$° 

+ £ a(<n 2 \(I a +U~ ^) ( {a) H 

a>l L 



A 



rak 
k 



(«) H ,m\ d 



The mixing for a = can be removed if we choose p to satisfy 

dt a(l +lna + p) = 



(72) 



(73) 



For a > 1, there are two conditions that must hold for the mixing to disappear, namely, 



I a + fa 



which can be easily solved to give 



fa 



A 



A + &m 2 a 
12 



A + 6m?, 



r 



(74) 
(75) 



(Note that the above expressions are free from singularities in view of (B16).) 

Thus, we see that the transformation (^) together with (E^) diagonalizcs the quadratic action as well as the 
boundary conditions, with /(r) determined from ( |73| ) and ([74|) . We note that, while a a , for a > 1, is determined 
from ( [75^ the coefficient an remains arbitrary. (This arbitrariness is related to the invariance of the theory under a 
gauge transformation of the graviton of the form flB15| ).) 

The diagonalized action now takes the form 



2M 3 r c / / dt (L(d) + L(H)) 



(76) 



where 



L(d) 



> Q>0 



M(t) „ 



a>l 



(77) 
(78) 



L(H) = I^ a ($«) 2 L PF (^H mk ;m 2 a ) 

. The properties of the diagonalized theory will be discussed 



a>0 



and 



and the functions K and M are defined in 
in more detail in the next section. 

Let us present here the final structure of the metric g m k{x,T) after all the transformations. Combining ([l6|), (|2^), 
(H), © and (H) we have 



g m k(x,T) = a(r) (l + d(x)pna +/(t)]) £ rofe (a;)+ 

+ (°) J ff mfc ( 2 ;)a> + 5(T)d ;mfe (x) + E (a) ^mfc(x)$ a (r) 



(79) 
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The functions /(r) and S(t) are given in terms of their coefficients (|73[), (74) and (|7J]) respectively. As we noted 
earlier, the function S(t) is determined up to an additive constant. The reason is quite clear. As we will discuss in 
appendix B, the massless field ^H mk is determined only up to a gauge transformation (B15) and, therefore, the 



term proportional to d ;m fe in (|54j) can be gauged away at any fixed point t = constant, in particular, on any one of 
the two branes, but cannot be gauged away in the whole bulk. 

IV. 4D EFFECTIVE THEORY 

To study further the properties as well as the physical implications of the diagonalized theory, we need to know the 
values of the coefficients present in ( |77| ) and ( f78| ) after integration over the fifth coordinate. Although so far we have 
treated both the cases, A = and A 7^ 0, on the same footing, it will be more convenient, in this section, to consider 
these two cases separately. 

A = 0: 



Let us note, at the outset, that A = does not automatically imply 9 mk — f]mk- F° r an arbitrary Einstein space, 
we have |TT| 

B, ma kb = -jg (irnblak ~ Irnkiab) + Qmakb ( 80 ) 

where C_ makb is the Weyl conformal tensor, and the solutions for any given A are completely degenerate with respect 
to different choices of C_ makb . 

When A = 0, it is clear from ( |74|) t hat f a = for a > 1, so that /(r) = p and the warp factor is given in pi]). 
With these, a direct evaluation of ([73[) determines 

2itkr r 

< (81) 



The coefficient in front of the kinetic term of the radion can now be determined from ( ]77j ) and (^54|) to be 

K d = dtK(r) = — / a(lna + p)(lna + p + 2)dt = -3irp > (82) 



2. 

Similarly, we can obtain the radion mass from d77|) and (|65|) and the explicit evaluation gives 



a 2 (\na + p)(\na + p+l)dt = ^2 a 2 — ^|(lna + p) 2 S(t - t t )dt = 

Therefore, the radion is massless, since 



_7r 1 = 0, 7T 



w MM , N 

— Y^dt = (83) 



The quadratic action, therefore, takes the form 

J d 4 x^l j -6Tr\p\M 3 r c d' m d, m - M 3 r c £ L PF ( {a) H mk -ml) J (84) 



Q>0 



Let us next normalize the fields d and ^ a 'H m k in the following way. Let us assume that exp(— 2nkr c ) <C 1 and express 
M in terms of the AD effective Planck mass Mj> L = M 3 /fc |J. Furthermore, let us define 



g fer c 



d(x) = —= D(x) 

V24irkr c M PL 

(B) flm*W = Tpirr {a) X mk {x) (85) 
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In terms of these fields, the quadratic action becomes: 



(86) 



The content of the theory, at the quadratic level, is now obvious - there is a massless scalar field (the radion), a 
massless spin-2 field (graviton) and an infinite tower of massive spin-2 fields with masses rr? a . We see that no ghost 
fields are present in this theory, since the kinetic term of the radion has the right sign, and the fields ^ a 'X mk carry 
precisely spin-2 content. We can, therefore, think of these fields as the physical fields. The non-presence of (ghost) 
scalar component in the fields ^X m k (for a > 1) is due to the conditions (B14) which these fields satisfy. Note that 
in our framework, the conditions (B14) arise naturally from the equations of motion for these fields (for more details 
see the appendix in |l9|]), and do not have to be imposed by hand. 

To understand the coupling of the fields D(x) and ^ a 'X m k to matter, we need to express the metric g m k(x, t) in 
terms of these physical fields. To that end, let us recall the expression for the normalized functions <& a (r) [|T| 



J2{x a ) 



s 2fcr c (r-7r) 



2 -2kr c 



'Y 2 (x a e 



kr c (r — TT) 



(87) 



for a/0, and $° — \Jkr c . Here, a: Q 's represent the positive zeroes of the Bessel function J\(x). Furthermore, the 
masses, m„, of the fields ^X m k are given by pit] 



2 _ j,2 2 — 2fcr e 7T 



On the visible brane, we have 



$"(T = 7T) 



, a > 1 



(89) 



while on the hidden brane they are rescaled roughly by a factor (a-dependent) of the order of exp(— 2kr c n). Collecting 
earlier results, we can now write the metric g m k(x, r) in terms of the physical fields (p3) as 



9mk{x,T) 



-2kr c T 



(y e -2fcr c 7T + T ) e fcr c 7r 



\/6itMpl 
1 



n>->))g mk (x)+ ^ T) fZ D -M*) 

-mk V247rfcr c AfpL 



Mpl \fkr c MpL 



(a) X mk {x)<$> a { T ) 



a>l 



(90) 



Of particular interest to us is the restriction of (|9(]) to the visible brane, which would determine the coupling of the 
fields D and ^ X mk to matter located there. 



9mk( x ) 



-2k 




(91) 
(92) 

(93) 



Some comments are in order. The couplings of the spin-2 states (|93|) were already discussed in [|2l|. The radion has 
two types of couplings to matter. The non-derivative one, on the visible brane, has strength comparable to that of the 
massive spin-2 states in (93) which is of the order of (TeV) -1 . This coupling is nonzero on the hidden brane as well, 
although much weaker (as also discussed in Iffl). The coupling on the hidden brane has its origin in the mixing of the 
radion with the graviton. The expression for the derivative coupling of the radion is nontrivial, due to the fact that 
we account for the presence of the massive KK tower of states. As we have noted earlier, the mixing of the radion, 
involving derivatives, with the graviton in (54) involves only the coefficient <to, which is a constant and can, therefore, 
be gauged away in the whole bulk. However, the mixing of the radion involving derivatives with the massive KK tower 
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in (pj) involves coefficients that depend on r, namely, the function S(t) in (|90j). As a result, in general, its value will 
be different at different points of the fifth dimension. The derivative coupling term can be arranged to vanish on the 
visible brane, but unless the function S(t) takes special values, the term D. m k cannot be simultaneously gauged away 
on both the branes. Let us estimate the coefficient of this term in (|92|). Since k < Mpl (see the number outside 
the parenthesis is of the order of (TeV)~ 3 . The term inside the parenthesis needs a more careful consideration. We 
use kr c — 12 and take the values for the first fifty roots, x a , with eight digit precision from A direct numerical 
evaluation shows that the first term of the series is —3.5593. The subsequent terms, in the series, are of alternating 
sign and tend to have very slowly decreasing magnitude. We believe that the series is absolutely convergent. The 
numerical evaluation of the series including the first fifty terms yields the value 0.74 times the first term. This value 
essentially stabilizes after approximately the twentieth term. The correction to this value from the remainder of the 
series is of the order of the next term, which is of the order of 1CP 3 . This is roughly the order of magnitude of each 
of the terms in the series for the next fifty or so terms. The 200-th term, in the series, for example, has a value of the 
order of 10 -4 . Therefore, we estimate the term in the parenthesis to be of the order of unity. 

A ^ 0: 



When A ^ 0, the diagonalized quadratic action is of the form 

S eff = J d*Xy/=l C(d) + J^x^lC 01 (94) 



a>0 



where 



and 



C(d) = Ci (d' m m ) 2 + C2 d> m d, m + Cm d 2 
C a = -M 3 r c L PF ( {a) H;m 2 a ) 



C m = -2 / dtM(r) (95) 

r c J --a 



( 2 - -2M 3 r, 

m 3 r 



The spin-2 content of the theory has properties similar to the case A = and does not need any further discussion. 
However, the radion now exhibits new features. The kinetic term has a fourth order derivative term whose coefficient 
has the sign of A. Although in principle the A ^ case can be analyzed in a fashion similar to the A = case, there 
are some technical complications which make the analysis quite difficult. For example, we do not know of a complete 



system of functions <f> Q that is of a simple form, although equation (B1C) can be solved in terms of associated Legendre 
functions or equivalently in terms of hypergeometric functions. Since the functions K(t) and M(t) (see (|64|), (pa)) 
are given in terms of the function f(r), which in turn can only be determined from its series expansion, we do not 
have much knowledge about these functions either. Therefore, we are unable to discuss the values of the coefficients 
( |95| ) any further. Nonetheless, it is easy to see that the coefficient £i is unlikely to be zero and has the same sign as 
A (It will be zero if all the I a, s are to vanish which is highly unlikely although not impossible.). This is sufficient to 
indicate that higher derivative ghosts will be associated with the radion field in this case. 



V. CONCLUSIONS 



In this paper, we have systematically studied, within the Lagrangian approach, the question of identification and 
consistent inclusion of the radion field in the two brane Randall-Sundrum model. We have exploited the symmetries 
of the theory to identify the radion in an unambiguous manner. Using the background field method, we have given 
a unified derivation of the classical solutions as well as the form of the action at quadratic order in the fields. We 
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have shown that the background metric, in general, corresponds to that of an Einstein space with a warp factor. We 
have discussed in detail the question of diagonalization of this quadratic action which is essential in understanding the 
spectrum of the theory. We have studied further the effective 4-dimensional action following from this diagonalized 
action and it is shown that the radion has no Kaluza-Klein tower and is massless for backgrounds with A = 0. In 
this case, the graviton as well as the infinite tower of states truly describe spin-2 particles and the theory is free of 
ghosts. For the case A ^ 0, however, the situation is less clear and higher derivative ghost terms, involving the radion, 
appear in the action. The question of matter coupling is also discussed. We would like to emphasize here that there 
are only two assumptions that we have made, i) the background r(x) — r c — constant and ii) a factorizable form of 
the background metric, g mn (x,T) — a(r)g (x). The entire analysis, otherwise, is quite general. 
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APPENDIX A: SOME USEFUL RELATIONS 



In this appendix, we collect some relations that are useful in the decomposition of the action as well as the junction 
conditions. For N defined in ([Toj) , and with (|l6|), (|l7|), we have the following decompositions up to second order in 
the fluctuations h mn and d. 

N = r 2 (l + 2d + d 2 -T 2 r 2 c g ab d, a d, b ) 



N 
1 

N 
1 



r c I I. + d + 'I 2 
1 



9 ao d, a d, 



l-2d + 3d 2 + t 2 t\ g ab d, a d :b ) 



1 



N 



l-d + d z 



g ah d,ad, b 



(Al) 



The determinant of the metric g mn can also be decomposed up to quadratic order as (see |15|) 

-h - -h a b h b „ + -I 
2 4 b a 8 



'-g = VS 1 



(A2) 



where h = g h m k, and we have, to quadratic order in the fluctuations, 



Vn \/—g = r c y/—g ( 1 



d+-h+ -dh 
2 2 



1 



1 



h\h b a 



c ~ab 



g ao d a d b 



(A3) 



The decomposition of RmL an d can be found in (note that the notations there differ slightly from ours). 



APPENDIX B: KALUZA-KLEIN DECOMPOSITION 



In this appendix, we describe the Kaluza-Klein decomposition for the fluctuation of the metric, H m k, and discuss 
some of the properties of the basis functions. Let us consider the action (f45|). When d = 0, it is easy to check from 
eqs. @ and @ that 



so that the quadratic action takes the form 

S(H) = -2A/ 3 r c J d 4 x^lj dt L(H) (Bl) 

where 
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L(H) 



a \ 



2 V 2 



-H a k;mH a 



tt Tjam:k i Tjmk tt 

tlak-.mtl +X1 



1 



-H m H, 



H' m k(H mk )'-(H 



4r 2 

oA 3 (a') 2 
8 + 1~rT 



E 



a 2 Vj 
r c 16M 3 



4 r : 

(2H mk H mk - H 2 ) 
(2H mk H mk - H 2 ) 5(t - U 



^ aa ' (o n' n mk u 1 tj\ 
~~ { 2H mk H -Ha.) 



Here, as usual, H = H g , and the field H m k has to satisfy the boundary condition (see (|53j) 



H mk {x,T)\ t=ti = 



Varying (Bl) with respect to H m k, we obtain 

H mk;a, a ~ H % 9_ mk + H- m k + H~ U ;ab g m) _ - i? ma;fc ' Q — H ka . m ° 



H m k - H g 



.ink 



Hmk — H g 



)--( 

mk J 4 \ 



2H mk -H g 



■ink 



= 



(B2) 
(B3) 
(B4) 
(B5) 

(B6) 



(B7) 



As we have seen from the discussion following (24), all singular terms from the equations of motion disappear when we 
use dot (r) derivatives. The above equation admits separation of variables and let us introduce a system of functions 
$ Q (r) such that: 



H mk (x,T) = {a) H mk {x)<S> a (T) 



(B8) 



Then, (|B7|) can be rewritten as 



(a) it ;a _ (a) TT-a , (a) tt , (a) rrab _ (a) tt ;a (a) it ;< 

n mk;a U ;o 9_ m k + • rnk + ; ab i-mk U ma-k n ka;m 



a 



(B9) 
(BIO) 



where represent the separation constants. 

Let us first analyze the system of functions <i> a (r), defined through equation ( BIO ). These functions satisfy the 
boundary condition 



<J> a (r) 







(Bll) 



as a consequence of the junction conditions ( |B6| ) and (B8). They define an orthonormal basis with respect to the 
scalar product: 



a{T)^ a {T)^{T)dt = S a[j 



(B12) 



The explicit forms for these functions, when A = 0, have already been obtained in || and (See also pi]| ). 

The equations (B^) are simply the equations of motion for spin-2 fields with mass m a on a background which is an 
Einstein-space JXS ] (for a detailed discussion of the flat case see also Q). In Q it is shown that, when m 2 a > 0, 



(a) 



H. 



2R 



'— ambk 
k 



il {a) H n 







{a) H = (a) H^ = 



(B13) 
(B14) 



and, therefore, we have five propagating degrees of freedom. Here, R am bk 1S the Riemann tensor constructed from the 
metric g mk - On the other hand, when m a = 0, the theory is invariant under the gauge transformation 



1G 



(B15) 



and describes a massless spin-2 field. In [Q it was also demonstrated that the masses m a of the states in the KK 
tower have to satisfy 



X + 6m 2 a 



(B16) 



Using (g8|) and (|B10|), (|34|), (|35|), ( |Bll| ) and (|B12|), we can rewrite the action (|Bl|) in terms of the KK states 
^Hjnk- After some algebra, we obtain 



S^(d = 0) = -2M 3 r c J2 J d^x^l J* dt |(<1> Q ) 2 L PF ( ^H mk ;m 2 a 



(B17) 



where Lpp is the generalization of the Pauli-Fierz Lagrangian for a background which is an Einstein space 
namely, 



Lpp{H,m 2 ) — -H a k- m H ak,m — H a k- m H 



am:fc i rrmfc 



1 



Fran tt TP-' m U 

■fe-n,m — n ., 



A 



(H m kH mk — H 2 ) + — (2H m i~H 



rak 



H 2 



(B18) 



It is easy to check that the equations of motion following from (B18) are precisely the same as (B9). 
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